Asymmetry of Nonlinear Transport and Electron Interactions in Quantum Dots 
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The symmetry properties of transport beyond the linear regime in chaotic quantum dots are 
investigated experimentally. A component of differential conductance that is antisymmetric in 
both applied source-drain bias V and magnetic field B, absent in linear transport, is found to 
exhibit mesoscopic fluctuations around a zero average. Typical values of this component allow a 
measurement of the electron interaction strength. 



Quantum transport in disordered mesoscopic conduc- 
tors and chaotic quantum dots has been widely stud- 
ied in the regime of linear response, and is well under- 
stood in terms of universal statistical theories pj, even 
in the presence of significant electron interaction 0, 0] ■ 
A central principle of linear mesoscopic transport con- 
cerns the symmetry of magnetoconductance: as a con- 
sequence of time reversal symmetry and microscopic re- 
versibility close to equilibrium, the differential conduc- 
tance g — dl/dV of a two-terminal sample is symmet- 
ric in magnetic field B: g(B) = g(-B) HIE 0, with 
generalized reciprocity (Landauer-Biittiker) relations for 
multi-terminal coherent conductors |fj . Beyond linear re- 
sponse, i.e., at sufficient applied bias that the current I is 
no longer proportional to source-drain voltage V, these 
symmetry relations break down. Unless disallowed by 
some special symmetry, differential conductance beyond 
linear response can generally contain a component (here 
denoted g) that is proportional to both B and V: 



g = aVB. 



(1) 



For instance, contributions to g of the form in Eq. 1 are 
permitted in non-centrosymmetric materials [jj , chiral 
conductors Q , including carbon nanotubes Q, or con- 
ductors with crossed electric and magnetic fields 10]. 

In the absence of electron interaction, including indi- 
rect interaction such as inelastic phonon scattering, the 
coefficient a in Eq. 1 vanishes, since conduction at each 
energy within the finite bias window independently obeys 
the symmetry of linear response and different energies do 
not mix 0,^2- Moreover, as discussed recently in Refs. 
[lii. fll ] , the coefficient a is proportional to the electron 
interaction strength. This suggests that terms of the form 
of Eq.^can be used to measure interaction strength. As 
pointed out in Refs. [l3l IT^|. finite g arises from quan- 
tum interference and is suppressed by decoherence and 
thermal averaging. 

This Letter presents a detailed study of the symme- 
try of nonlinear conductance in an open chaotic GaAs 
quantum dot, using gate-controlled shape distortion to 
gather ensemble statistics. We focus particular attention 
on the component of g that is odd in both B and V, 
as in Eq. ^ at moderately small source-drain voltage, 



V < A/e, and fields, B < <Po/A, applied perpendicular 
to the plane of the dot, where A is the average quantum 
level spacing in the dot, 0o = h/e is the flux quantum, 
and A is the dot area. We find that the component of dif- 
ferential conductance antisymmetric in B, denoted gs-, 
is also largely antisymmetric in V and shows mesoscopic 
fluctuations as a function of B, V, and shape-defining 
gate voltage Vq. As anticipated theoretically |13Lll4j . we 
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FIG. 1: (a) Differential conductance g for N = 2 as a func- 
tion of magnetic field B (solid curves) and —B (dashed curves) 
at source-drain bias V as indicated, at base electron temperature 
T e i = 45 mK. The blue (red) curves are offset by +0.5 (+1 ) e 2 /h, 
respectively, (b) Antisymmetric in B conductance gB- of the 
traces shown in (a), as a function of B at V (solid curves) and 
— V (dashed curves) as indicated. The blue (red) curves are off- 
set by +0.1 (+0.2) e 2 /h, respectively. Black curves demonstrate 
g(B,V = 0) = g(—B, V = 0) and blue/red curves indicate gs— is 
largely antisymmetric in V. 
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find that the average coefficient a measured over an en- 
semble of dot shapes vanishes. The standard deviation of 
a, denoted 5a, does not vanish and is used to character- 
ize the strength of the interactions, as discussed below. 
The dependence of 8a on the number of modes N in 
the quantum-point-contact leads is found to be in dis- 
agreement with theory However, present theory 
assumes that electrons do not thermalize within the dot 
and that decoherence effects are negligible. Accounting 
for an increased amount of thermalization within the dot 
at smaller N (consistent with independent measurements 
of electron distribution functions 15]) appears able to ex- 
plain qualitatively the observed dependence of da on N. 

Measurements were carried out using a quantum dot 
of area A ~ 1 /an 2 formed by Ti/Au depletion gates 
[see Fig. EJd)] on the surface of a GaAs/Al .3Gao.7As 
heterostructure 105 nm above the 2D electron gas. A 
bulk electron density n ~ 2 x 10 11 cm~ 2 and mobility 
/i~2x 10 5 cm 2 /Vs, giving a mean free path £ ~ 1.5 jum, 
indicates ballistic transport within the dot. This device 
contains Ndot ~ 2000 electrons and has an average level 
spacing A = 2irh 2 /m*A ~ 7/^eV, where m* — 0.067 m e 
is the effective electron mass. The dot was designed to 
lack spatial symmetry [l6j and is found to show univer- 
sal statistics in linear conduction characteristic of chaotic 
classical dynamics 0, 0, El ■ The gate marked "shape" , 
with voltage Vg applied, was used to modify the dot 
boundary , p roviding a means of generating ensemble 
statistics |l7[ . Gates n, vol, and w2 were actively modu- 
lated depending on Vg to maintain an integer number of 
quantum modes N in each lead throughout the ensem- 
ble of shapes. Gates pi and p2 were strongly depleted, 
isolating the adjacent smaller dot, which was used only 
to measure the electron temperature (T 6 ;= 45 ± 5 mK at 
base) using Coulomb blockade peak width. 

To measure g, simultaneous lock-in measurements of 
both differential current and voltage across the dot were 
made using a four-wire setup at 94.7 Hz with typical 
ac excitation of 2 fj,V, in the presence of a dc source- 
drain voltage V. All voltages were applied symmetrically 
across the dot [Fig. EI a)] to reduce self-gating [lfl |. The 
measurement sequence was as follows: g(B,V, Vq) was 
measured as a function of V (innermost loop) and B (2nd 
loop, with higher point-density around zero field). Each 
of these two-parameter sweeps was repeated 20 times to 
reduce noise (3rd loop), for each of 16 statistically inde- 
pendent shape gate voltages Vg (4th loop). Each multi- 
parameter sweep, which took about 30 h to complete, was 
measured at N = 1,2, 4 modes in each lead (outermost 
loop). 

As expected for a two terminal device 0, 0, @], the 
linear conductance g(B, V = 0) was found to be sym- 
metric in B within measurement resolution, as seen in 
Fig-^a). At finite bias, however, the ±B symmetry is 
broken, as seen in the blue and red traces in Fig. [Ja). 
Decomposing g into components that are symmetric 
((Jb+) and antisymmetric (gB-) in field, gB±{B,V) = 
\g{B,V)±g(-B,V)) /2, we find g B -(B,V = 0) Goin- 
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FIG. 2: (a) Antisymmetric component gs— (color scale) as a func- 
tion of B and V at N = 2. Dashed lines indicate cuts shown in (b) 
and (c). (b) gB— as a function of B at V = 0, 7.5, 15, 55 fiV. (c) 
gB— as a function of V at B = 0, 0.3, 0.8 mT. Linear fits to curves 
in (b) and (c) are shown in dashed blue lines, with solid blue lines 
indicating the fitting range. Inset to (c) shows the same cuts as 
in (c) with the V axis extended to ±125 fiV and the gB- axis ex- 
tended to ±25 X 10~ 3 e 2 /h. (d) Electron micrograph of a device 
with the same gate design as the one measured. Only the yellow 
gates are used in this experiment. V is applied symmetrically. 



parable to the measurement noise of 3 x 10~ 3 e 2 /h and 
barely visible, while at V = ±20/ieV and ±200^eV, 
gs- is sizable and is found to be largely antisymmet- 
ric in V. The typical amplitude of fluctuations of gs- is 
~ 0.025 e 2 /h, much smaller than the amplitude of fluc- 
tuations of g, which is ~ 0.5e 2 /h. 

The antisymmetry of gs-(B, V) in source-drain volt- 
age V is also seen in Fig. 0{a), which shows gB-(B,V) 
for N — 2 as a two dimensional color plot for a differ- 
ent shape sample than the one shown in Fig. ^ Note 
that gB-(B, V) changes sign repeatedly as a function of 
both B and V . Similar characteristics are seen for all 
dot shapes and values of TV, though for some shape sam- 
ples the antisymmetry in V is more pronounced than in 
others. 

For small fields, B < <fio/A, gB-(B) is proportional 
to B and becomes sizable for V > A/e ~ 7 fieV, as 
seen in Fig. Hfb). At larger fields, B > <f> /A, g B -(B) 
shows mesoscopic fluctuations as a function of B, includ- 
ing sign changes. The field scale where this crossover oc- 
curs is consistent with the field scale of weak localization 
[see Fig.|3Ja)] and conductance fluctuations in this device 
(not shown). This field scale is somewhat smaller than 
(j>o/A ~ 4mT, similar to previous experiments |l8lll9| be- 
cause the relevant area is not the dot area but the area of 
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typical trajectories before escape, which is larger than A 
by a factor ~ </A^ HJ. For V < A/e, g B -{V) is also 
found to be proportional to V and becomes sizable for 
B > lmT, as seen in Fig.^c). For V > A/e, g B -{V) 
starts to deviate from the linear dependence on V and 
shows mesoscopic fluctuations as a function of V, also 
including sign changes. As mentioned above, while some 
shapes show more pronounced antisymmetry of gs- with 
respect to V than others, g B - is predominantly antisym- 
metric in V, see discussion of Fig.QJa). The characteris- 
tic scales V ~ A/e and B ~ <fio /A for mesoscopic fluctu- 
ations of g B -(B,V) provide assurance that this compo- 
nent arises from coherent transport within the dot. 

The average symmetric component of conductance, 
(dB+)v G i which is maximal around V = 0, has a pro- 
nounced minimum around B = due to ballistic weak 
localization [2lj], as seen in Fig. EH a). From the magni- 
tude of the weak localization feature in (g(B, V — 0))y G , 
we extract a base-temperature phase coherence time 
t v ~ 2 ns, consistent with previous experiments |2l|. As 
seen in Fig.^a), finite V reduces the dip (g(B, V = 0))y G 
around B = 0, presumably the result of dephasing caused 
by heating and non-equilibrium effects [22L |23| . Fig- 
ure Efb) shows that the standard deviation of the sym- 
metric component 5g B+ is peaked around B = 0, a 
well-known effect associated with the breaking of time- 
reversal symmetry 0, 0, 0, • Similar to the weak lo- 
calization correction, Sg B+ is also reduced at finite bias. 
In this case, the reduction is due both to dephasing and 
an explicit dependence on temperature and energy aver- 
aging. 

Figure 3(c) shows the shape-averaged antisymmetric 
component of conductance {g B -)v G , which is zero within 
statistical uncertainty (based on 16 shape samples), in- 
dicating that fluctuations of g B _ cancel each other upon 
averaging, as predicted theoretically [l3L IT^ | . As seen in 
Fig- Eld), the standard deviation Sg B - becomes sizable 
on a voltage scale of order of the dot level spacing, with a 
maximum at moderate |V|, and then decreases for larger 
\V\, presumably due to heating or decoherence effects. 
By definition, 5g B - is zero for B = and, as expected 
for a g B - that is largely antisymmetric in V, Sg B - is 
comparable to the statistical error along the V = line. 

To characterize symmetry in terms of statistical 
quantities, we define normalized symmetry parameters 
C B ±v± = 6 2 g B ± v ±/5 2 g B ±, where S 2 g B ± v ± is the vari- 
ance of the component of g with a particular symme- 
try with respect to B and V. Here, ensemble averaging 
is performed over the entire (B, V, Vg) parameter space, 
resulting in approximately 5 x 4 x 16 statistically inde- 
pendent samples. We note that C ranges from zero to 
one, and C B ±v± + C B ±y T = 1. As seen in Fig. 0] 
C B +v+ ~ 1 (open squares) and C B +v- ~ (filled 
squares) within the statistical uncertainty, showing that 
g B + is symmetric in V for all measured conductance sam- 
ples and modes numbers N. C B -v- ~ 0.75 (open cir- 
cles) and C B -v+ ~ 0.25 (filled circles) without a strong 
N dependence, indicating that a significant part of g B - 
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FIG. 3: Average conductance component (gB+)vc ( a ) an( ^ 
{9b-)v g ( c ) an d standard deviation Sgg + (b) and Sgg_ (d) (color 
scale) obtained from 16 independent shape samples at N = 2 modes 
as a function of B and V. (gB—) is zero within statistical uncer- 
tainty; standard deviations are used to characterize the antisym- 
metries. (gB+) shows a strong dip at (B, V) = due to ballistic 
weak localization. 



is antisymmetric in V, though the antisymmetry is not 
perfect. 

The interaction coefficient a was extracted for each 
shape by performing fits linear in both B and V to 
g B _(B,V) for small V < A/e and B < (j>o/A, similar 
to the fits shown in Figs. Efb) and (c). The resulting a 
is seen to have mesoscopic fluctuations and changes sign 
frequently, consistent with (a)v G — within statistical 
uncertainty. The standard deviation da, used to char- 
acterize the strength of electron interaction, is shown in 
Fig. IHb) (open triangles) as a function of the mode num- 
ber TV. 

For a comparison of the interaction parameter a to 
theory, we note that Ref. [14| predicts 5a for kT <C A^A, 
V < A/e and B < <j> Q /A: 



8a = 8a' 



1 



A e 2 



2A 2 A. 



(2) 



where 6a' is a dimensionless parameter characterizing the 
electron interaction strength. Using a different approach, 
Ref. |13| arrives at a similar expression for 5a, with the 
same dependencies on N and A. To facilitate a compar- 
ison of experiment and theory, we plot the dimensionless 
quantity 5a' in Fig. Elb), noting a pronounced depen- 
dence of 5a' on N. As defined in Eq. 2, 5a' is not ex- 
pected to depend on N, according to Refs. [Tol Il4|. A 
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FIG. 4: (a) Normalized symmetry parameter C as a function of 
the number of modes N for all B and V symmetries as indicated, 
giving a statistical characterization of the nonlinear field symme- 
try properties of conductance over all V , B and shape samples 
measured, (b) Standard deviation 8a (open triangles) and dimcn- 
sionlcss standard deviation 8a' (solid diamonds) of the interaction 
parameter (see text) as a function of TV. 



likely reconciliation is that for N < 4, the escape time 
from the dot is sufficiently long so that electrons have 
time to equilibrate (thermalize at a higher temperature) 
before escaping. As a result, 6a' is reduced at N — 1, 2 
compared to N = 4. Indeed, from direct measurement 
of electron distribution functions in the same device |15| , 
it is known that for N — 1 the distribution of electron 
energies in the dot is thermal, with an elevated effective 
temperature that depends on V. For N = 4, on the other 
hand, with a shorter electron dwell time in the dot, non- 
thermal distributions are seen ^{|. For the case N = 4, 
the measured 6a' — 1.1 ± 0.2 is reasonably consistent 



with the theoretical estimate 6a' — n ^| , which assumes 
perfect screening of electrons. Other factors possibly con- 
tributing to the observed N dependence of 6a' may in- 
clude imperfectly transmitting modes or mode mixing as 
well as finite temperature and/or decoherence effects, all 
of which are not currently accounted for theoretically. 

In conclusion, we have investigated the magnetic field 
asymmetry of conductance beyond the linear regime in 
gate-defined quantum dots. The conductance component 
gs-, which is defined to be antisymmetric in B, was 
found to be also predominantly antisymmetric in V and 
is of the form of Eq. [T] for V < A/e and B < <t>o/A. 
The interaction coefficient a, extracted from linear fits 
to gs-, has mesoscopic fluctuations with zero average. 
Comparison to recent theory 0, ^| is most appropriate 
for the data with four modes per lead, where electrons 
remain out of equilibrium during the short dwell time. In 
this case, consistency with theory 0] relating a to the 
interaction strength appears reasonable, suggesting good 
electron screening. 
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Note added. During completion of this manuscript, 
closely related experimental work on carbon nanotubes 
appeared |24|. 
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